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We inaugurate a framework for studying preheating and parametric resonance after inflation with- 
out resorting to any approximations, either in gravitational perturbation theory or in the classical 
evolution of the field(s). We do this by numerically solving the Einstein field equations in the 
post-inflationary universe. In this paper we show how to compare our results to those of gauge in - 
variant perturbation theory. We then verify Finelli and Brandenberger's analysis (hep-ph/9809490) 
of super-horizon modes in m 2 2 inflation, showing that they are not amplified by resonant effects. 
Lastly, we make a preliminary survey of the nonlinear couplings between modes, which will be 
important in models where the primordial metric perturbations undergo parametric amplification. 
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I. Introduction Postulating a period of nearly expo- 
nential growth in the primordial universe, inflationary 
cosmology JIJ solved many problems which plagued pre- 
vious models of the big bang. Immediately after inflation 
the universe is dominated by the inflaton, the scalar field 
whose evolution controls the dynamics of the inflationary 
era. The inflaton energy density is then converted into 
thermalized, "normal" matter, a process called reheating. 
The first predictions were that the universe would 
be reheated to a temperature several orders of magnitude 
(or more) below the inflationary scale. Subsequently, it 
has been realized that non-equilibrium, non-perturbative 
processes can drive explosive particle production at the 
end of inflation [p|df9| . This initial era, dubbed preheat- 
ing, is governed by Mathieu-like equations, and the rapid 
particle creation associated with the instability bands of 
these equations is often termed parametric resonance. 

Developing a full understanding of preheating has been 
a long and arduous task. There were two main compli- 
cations. Firstly the difference between parametric reso- 
nance in Minkowski space and in an expanding universe 
is profound. In the latter case, the instability bands may 
be time-dependent and so a particular mode may pass in 
and out of resonance. Secondly, the back-reaction of the 
created particles onto the other fields must become sig- 
nificant; and indeed the parametric resonance will even- 
tually terminate. With a few exceptions though, (see 
e.g. ]20|-|23|]), spacetime was described by the usual Fried- 
man Robertson Walker metric. Although the scale fac- 
tor could be determined self-consistently (see e.g. [p^[), 
the possibility of resonant growth of metric perturba- 
tions was not considered. However, in a recent paper 
Bassett et al. |22| highlight the importance of couplings 
between perturbations in the fields and in the metric. As 
they point out, metric perturbations are sources for the 
field perturbations, and vice versa, so it cannot be as- 
sumed that the role of metricperturbations in preheat- 
ing is only slight (see also |2p| , |2l| , ^3|] ) . During preheating, 
in this framework, parametric resonance will necessarily 
amplify metric perturbations and so nonlinear gravita- 
tional effects will become significant. 

In this paper, we abandon the perturbative approach 
entirely and study the post-inflationary era by numeri- 



cally solving the inhomogeneous, nonlinear Einstein field 
equations, and derive the connection between our results 
and the predictions of perturbation theory. We then con- 
firm the conclusions of Finelli and Brandenberger |23| ] 
who use gauge invariant perturbation theory to show that 
super-horizon modes are not parametrically amplified af- 
ter m 2 <f> 2 inflation. Lastly, we use the m 2 </> 2 theory as a 
toy model to explore the consequences of having signifi- 
cantly amplified modes in the post-inflationary universe. 
In particular, we see how power can be transferred from 
modes with large amplitudes to all other modes. The 
mode-mode coupling thus appears to broaden the impact 
of resonance in the early universe. This is the first time 
in the study of preheating that the mode-mode coupling 
of the field and metric perturbations has been calculated 
in a fully relativistic manner. In m 2 <p 2 inflation, met- 
ric perturbations are not amplified so we inject a large 
perturbation via our initial conditions. In future work 
we will give a nonlinear treatment of inflationary models 
such as Xcj) 4 and g 2 <fi 2 x 2 , where the field perturbations 
are amplified by resonant effects, and study the role of 
mode-mode coupling in these systems. 

II. Metric and Field Equations We consider scenarios 
whose initial state corresponds to a spatially flat Fried- 
man Robertson Walker universe with a perturbation con- 
sisting of a single Fourier mode. The wave-vector, k, of 
this mode can be assumed to be parallel to the z direc- 
tion, so we adopt a metric where the x and y directions 
are homogeneous and the inhomogeneous contributions 
are functions of z alone. The simplest choice of metric 
which expresses this constraint is 

ds 2 = alt 2 - A 2 (t, z) dz 2 - B 2 (t, z) (dx 2 + dy 2 ). (1) 

A single scalar field has the Lagrangian 
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and equation of motion 
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Overdots indicate derivatives with respect to t, and 
primes with respect to z. The non-zero components of 
the Einstein field equations, G^ v — — k 2 T m „, are therefore 
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Equations (||) and (Q) are constraints which must be sat- 
isfied by the initial conditions, and are then used to check 
the accuracy of our numerical integration. 

We choose A, B and to be independent of z at t = 0, 
so the initial inhomogeneity is expressed by A, B and 0. 
This makes it easier to ensure that the initial data satis- 
fies the constraints, but does not involve an unacceptable 
loss of generality. With A'(0,z) = B'(0, z) = 0'(O, z) = 
the first constraint, equation (g), is trivial. Rescaling the 
spatial co-ordinates fixes ^4(0, z\ = B(0, z) = 1, while we 
denote 0(0, z) — 0q. Equation (m) requires 
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where we pick C so that (A(0,z)) — (B(Q,z)), where 
(■ • •} is a spatial average. This normalizes the metric 
functions and scales the co-ordinates, while 0(0, z) deter- 
mines the configuration of the primordial universe when 
we begin our simulation. All simulations presented here 
begin at the moment when inflation ends in the unper- 
turbed model, defined by d(t) = 0, where a(t) is the scale 
factor of the unperturbed universe. 

To facilitate comparison with the perturbative result, 
we excite a single /c-mode 
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where Z is the length of our "box" . This definition 
rescales A; in units of Z, and we use this form throughout 
the rest of this paper. The magnitude of the perturba- 
tion is fixed by e. Most inflationary models predict that 
eCl, but we can consider an arbitrary e, since we are 
working with the full Einstein equations. 

We use solve the field equations on an N point spatial 
grid. The simulations presented here all have N = 1024. 
We assume periodic boundary conditions, and use a 
fourth-order differencing scheme to express equations (Q) 
and (|5|) as a set of N ordinary differential equations, 
which are solved with a fourth-order Runge-Kutta in- 
tegrator. The resulting system suffers from a buildup 
of short- wavelength numerical noise, common in free- 
evolution schemes such as ours, due to the discretization 
process which can only approximate the true time-like 
hypersurface from timestep to timestep |24j] . To control 
this we use a pseudo-spectral approach, de-aliasing [ p5[ . 
This is implemented by filtering A, B, (and their as- 
sociated velocities), removing Fourier components with 
wavelengths less than 8 grid spacings at each timestep. 
Thus our effective spatial resolution is ~ 16 times larger 
than the separation of the individual points, Az. 

III. Testing Perturbation Theory Our first goal is to 
test the perturbative analysis of the interaction between 
the field and background perturbations in an inflation- 
ary model with the potential, V(<f>) = m 2 2 /2 [|o)-|3). 
In particular, Finelli and Brandenberger showed that 
first order gravitational perturbation theory does not pre- 
dict the parametric amplification of scalar perturbations 
with wavelengths considerably larger than the Hubble 
radius, as was suggested by Bassett et al. [p2[ . 

To independently check Finelli and Brandenberger's 
prediction and to confirm that perturbation theory gives 
an accurate description of the universe at the end of in- 
flation, we use our nonlinear code to analyze the same 
system. Consequently, we must relate our metric, which 
corresponds to a particular gauge choice, to gauge invari- 
ant cosmological perturbations. 

Working in physical time, with our symmetry con- 
straint and choice of flat spatial sections in the unper- 
turbed universe, the general scalar perturbation is p6| , 

ds 2 = (1 + 2tp(t, z))dt 2 - 2a(t)b(t, z), i dx i dt 

~a 2 {t) [(1 - 2if>(t, z))Sij + 2E(t, z) tij ] dx'dx (11) 

where commas denote ordinary derivatives. The back- 
ground solution, a(t), is defined by the FRW universe 
with our chosen initial conditions and e = 0. Comparing 
with equation (hi), we deduce that 



<p = 0, b' = 
A 2 =a 2 (t){l-2iJj + 2E"), 
B 2 = a 2 {t)(l-2ijj). 
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We are now in a position to calculate the gauge invariant 
metric perturbation $, 
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FIG. 1. The evolution of the <E>i mode, with a wavelength 
Z, is plotted for e = 0.001, After an initial transient decays 
the solution tends towards a constant, as expected for a su- 
per-horizon mode. The time is normalized so that m — 1. For 
our choice of phase, only the imaginary part of $1 is excited. 
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where we have set integration constants in the final ver- 
sion to ensure that $ vanishes when there is no pertur- 
bation. We can express E" and, therefore, <I>" in terms 
of the metric functions and their derivatives 



$" = BB + B 2 - AA - A 2 + 

H(A 2 - B 2 ) + 2H(AA - BB), 



(16) 



where H = a /a. Perturbation theory predicts the Fourier 
modes, of <£>. Recalling the Fourier transform, J 7 , 
of a derivative, we obtain = — fc _2 JT($") fc . How- 
ever, in linear perturbation theory $k is obtained from 
Mukhanov's variable, Q p7p , which satisfies p6| , p0 \ 
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Qk = 0, (18) 
(19) 



We ran our code for an initial perturbation with a 
wavelength equal to the size of our simulation region, 
which is 10 2 larger than the post-inflationary Hubble 
length, and e = .001. This value of e is somewhat ar- 
bitrary, but corresponds to a density contrast of roughly 
the same magnitude as those observed at COBE scales. 
We extracted as described above, and solved numeri- 
cally for Q with initial conditions chosen to match those 
of the nonlinear simulation. We find that the nonlin- 
ear and the perturbative results agree to within approxi- 
mately 1 part in 10 5 . The evolution of $k corresponding 
to the single excited mode is shown in Fig 1. This verifies 
both the perturbative treatment of the post-inflationary 



FIG. 2. The evolution of |$s| (upper) and |$io| are plot- 
ted, for the nonlinear case of e = .1. In this case, there is a 
significant transfer of power to higher modes. 



era in the absence of resonance and the specific conclu- 
sions of Finelli and Brandenberger pq ], i.e. that for 
m 2 (f> 2 inflation, is constant for super-horizon modes 
(15) and does not undergo parametric amplification. 



IV. Beyond Linear Gravity We now consider the role 
of nonlinear gravitational effects during preheating. A 
full relativistic treatment (or higher order perturbation 
theory) couples Fourier modes which are independent at 
first order. Post-inflationary perturbations are small, but 
in 4 inflation (and in more complicated models), para- 
metric resonance can increase $ dramatically. If this 
growth lasts long enough, nonlinear gravitational effects 
are unavoidable. 

We leave the full numerical analysis of models with 
more complicated potentials for future work, and use our 
m 2 (f> 2 system as a toy model to investigate the conse- 
quences of one or more modes acquiring a large ampli- 
tude. Since we have just shown that super-horizon modes 
are not amplified with a 4> 2 potential, we must inject this 
through our initial conditions. 

We examine a perturbation with k = 5, and e = 0.1, 
which is large in contrast to the primordial COBE scale 
perturbations, but is not inconsistent with a mode am- 
plified by parametric resonance. Two are plotted in 
Fig. 2. Modes for which k is an integer multiple of the k 
of the initial perturbation are excited. The k = 5 mode 
is outside the initial Hubble radius, but nonlinear effects 
cause it to grow, and the perturbative prediction is not 
applicable. We also observe that in initially over-dense 
(with respect to a flat, unperturbed universe) regions, A 
begins to decrease, signifying the onset of gravitational 
collapse. 

Finally, we consider a more complex case where two 
modes, whose k values are not multiples of one an- 
other, are initially excited. Specifically, we replace equa- 
tion©by 
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FIG. 3. Several modes of $ are plotted with e = 0.05 and 
k = 8, and the perturbation given by equation (pp|). Here, 
$i is initially excited, but significant amounts of power are 
transferred to other modes. 



and unamplified by any resonant processes. We also 
made a preliminary exploration of interactions between 
modes when the perturbation is large, demonstrating the 
ability of nonlinear effects to couple modes, which allows 
a transfer of power from modes which are parametrically 
enhanced to those which are not in resonance at the per- 
turbative level. Consequently, we conjecture that the 
nonlinear gravitational effects will broaden any resonance 
bands seen in the conventional perturbative analysis; we 
are currently investigating these questions in detail. 
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The two modes we initially excite become mixed in the 
gauge invariant perturbation, so $i always has a signifi- 
cant amplitude. This can be seen in Fig. 3, which corre- 
sponds to the choice e = 0.05 and k = 8. However, unlike 
the previous case where only modes which were multiples 
of the initially excited mode were enhanced, we now see 
a transfer of power between the initially excited "band" 
and a wide variety of other modes. Consequently, if some 
modes acquire a significant amplitude through paramet- 
ric resonance, nonlinear effects can then transfer power 
to modes which fall outside the resonance bands. 

V. Discussion Nonlinear effects may be important in a 
number of ways not considered here. For instance, when 
inflation is driven by two interacting scalar fields, the ho- 
mogeneous post-inflationary universe may exhibit chaos 
p8[ , and the interaction between this chaos and nonlinear 
inhomogeneous effects is intriguing. Likewise, paramet- 
rically enhanced perturbations will grow further due to 
nonlinear gravitational effects and could eventually un- 
dergo gravitational collapse, leading to the formation of 
primordial black holes in models where their existence 
is not currently predicted J2^|. The scattering of scalar 
perturbations is also a nonlinear effect and will lead to 
gravitational wave production 1 29 3(],|2^] . 

By abandoning perturbation theory and solving the 
Einstein field equations directly, we have laid the foun- 
dation for a fully nonlinear analysis of preheating in an 
expanding universe, which will include all possible cou- 
plings between perturbations in spacetime and the classi- 
cal field(s). We verified Finelli and Brandenberger's anal- 
ysis of m 2 4> 2 inflation, showing that the super- horizon 
modes of the gauge invariant perturbation are constant, 



[1] A. H. Guth, Phys. Rev. D 23, 347 (1981). 

[2] L. F. Abbott, E. Farhi, and M. B. Wise, Phys. Lett. B 

117, 29 (1982). 
[3] A. Albrecht, P. J. Steinhardt, M. S. Turner, and F. 

Wilczek, Phys. Rev. Lett. 48, 1437 (1982). 
[4] A. D. Dolgov and A. D. Linde, Phys. Lett. B 116, 329 

(1982). 

[5] J. H. Traschen and R. H. Brandenberger, Phys. Rev. D 

42, 2491 (1990). 
[6] L. A. Kofman, A. D. Linde, and A. A. Starobinsky, Phys. 

Rev. Lett. 73, 3195 (1994). 
[7] Y. Shtanov, J. Traschen, and R. Brandenberger, Phys. 

Rev. D 51, 5438 (1995). 
[8] D. Kaiser, Phys. Rev. D 53, 1776 (1996). 
[9] S. Y. Khlebnikov and I. I. Tkachev, Phys. Rev. Lett. 77, 

219 (1996). 

[10] E. W. Kolb, A. Linde and A. Riotto, Phys. Rev. Lett. 

77, 4290 (1996). 
[11] S. Y. Khlebnikov and I. I. Tkachev, Phys. Rev. Lett. 79, 

1607 (1996). 

[12] S. Y. Khlebnikov and I. I. Tkachev, Phys. Lett. B 390, 
80 (1996). 

[13] D. Boyanovsky, H. J. de Vega, R. Holman, and J. F. J. 

Salgado, Phys. Rev, D 54, 7570 (1996). 
[14] T. Prokopec and T. G. Roos, Phys. Rev. D 55, 3768 

(1997). 

[15] D. Boyanovsky, D. Cormier, H. J. de Vega, R. Holman, A. 
Singh, and M. Srednicki, Phys. Rev, D 56, 1939 (1997). 

[16] S. A. Ramsey and B. L. Hu, Phys. Rev. D 56, 678 (1997). 

[17] L. Kofman, A. Linde, and A. Starobinsky, Phys. Rev. D 
56, 3258 (1997). 

[18] P. B. Greene, L. Kofman, A. Linde, and A. A. Starobin- 
sky, Phys. Rev. D 56, 6175 (1997). 

[19] B. R. Greene, T. Prokopec, and T. Roos, Phys. Rev. D 
56, 6484 (1997). 



4 



[20] H. Kodama and T. Hamazaki, Prog. Theor. Phys. 96, 
949 (1996). 

[21] Y. Nambu and A. Taruya, Prog. Theor. Phys 97, 83 

(1996) . 

[22] B. Bassett, D. Kaiser, and R. Maartens, friep-ph/9808404 
(1998). 

[23] F. Finelli and R. Brandenberger, [hep-ph/980949C| (1998). 

[24] M. W. Choptuik, Phys. Rev. D 44, 3124 (1991). 

[25] D. Gottlieb and S. Orszag, Numerical analysis of spectral 

methods: theory and applications (SIAM, Philadelphia, 

1977). 

[26] V. Mukhanov, H. Feldman, and R. H. Brandenberger, 

Phys. Rep. 215, 203 (1992). 
[27] V. F. Mukhanov, Sov. Phy s. JETP 67, 12 97 (1988). 
[28] R. Easther and K. Maeda, |gr-qc/97110"35| (1997). 
[29] S. Khlebnikov and I. Tkachev, Phys. Rev. D 56, 653 

(1997) . 

[30] B. A. Bassett, Phys. Rev. D 56, 3429 (1997). 



5 



